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Recent ARPES measurements on the electron-doped cuprate Nd2_iCei,Cu04 can 
be interpreted in a mean field model of uniform doping of an antiferromagnet, with 
the Mott gap closing near optimal doping. Mode coupling calculations confirm the 
mean field results, while clarifying the relation between the Mott gap and short- 
, range magnetic order. The same calculations find that hole doped cuprates should 

■ follow a strikingly different doping dependence, involving instability toward spiral 

Q ' phases or stripes. Nevertheless, the magnetic order (now associated with stripes) 

^ ' again collapses near optimal doping. 



1. Introduction 
\Q ' 1.1. Mode Coupling Theories 

QQ , In conventional itinerant ferro- and antiferromagnets, mode-coupling 

' theories^^^ have proven of value in treating the role of fluctuations in reduc- 

. ing or eliminating long-range order, as well as in the development of local 

-« \ moments. Curie-like susceptibility, and in general the crossover to mag- 

' netic insulators. Similar approaches have been applied to charge-density 

wave systems'^ and the glass transition^. 

Attempts to apply such a formalism to study the antiferromagnetism 
of the cuprate superconducting compounds have been frustrated, since the 
Q \ antiferromagnetic (AF) phase is found to be unstable against hole doping, 

O ' toward either an incommensurate AF phase^ or phase separation^. Here, it 

is demonstrated that this tendency to instability is either absent or greatly 
reduced in electron-doped cuprates, and the mode coupling analysis can 
;J] ' provide a detailed description of the collapse of the Mott gap with doping. 

The results are of great interest of themselves: the collapse is associated 
with one or more quantum critical points (QCPs), and superconductivity 
is optimized close to one QCP. However, the results have an additional 
importance in the light they shed on the more complicated problem of 
the hole doped cuprates. First, the effective Hubbard U parameter has a 
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significant doping dependence, from U ^ W at half filling (where W = 8t is 
the bandwidth and t the nearest neighbor hopping parameter) to ?7 ~ W/2 
near optimal doping. Secondly, the same mode coupling theory which works 
so well for electron-doped cupratcs, breaks down for hole doping, due to 
the above noted electronic instability. This suggests that (a) the electron- 
hole asymmetry must be due to a band structure effect and (b) a suitable 
generalization of mode coupling theory should be able to incorporate the 
effects of this instability. 



1.2. ARPES of Nd2-xCe^CuOA 

Recently, Armitage, et al.^ succeeded in measuring angle-resolved photoe- 
mission spectra (ARPES) of Nd2-i;Cca;Cu04 (NCCO) as a function of 
electron doping, from essentially the undoped insulator x = to opti- 
mal doping X = —0.15. The doping dependence is strikingly different from 
that found in hole-doped La2- a;Sra;Cu04 (LSCO)*^'^. Both systems start 
from a Mott insulator at half filling, with ARPES being sensitive only to 
the lower Hubbard band (LHB), approximately leV below the Fermi level 
Ep- With hole doping, the LHB remains well below Ep while holes are 
added in mid-gap states (as expected, e.g., in the presence of nanoscale 
phase separation^°); for electron doping, the Fermi level shifts to the up- 
per Hubbard band (UHB), and the electrons appear to uniformly dope the 
antiferromagnet . 

Remarkably, the full doping dependence can be simply described by a 
mean- field (MF) t — t' — U Hubbard model^^, where t' is the second neighbor 
hopping, U is the onsite coulomb repulsion, and a one band (copper only) 
model was assumed for simplicity." A key finding is that the Hubbard 
U is doping dependent, leading to a quantum critical point (QCP) just 
beyond optimal doping, where the Mott gap vanishes. The appearance 
of a peak in superconductivity near an AFM QCP is a fairly common 
occurance'^^; in particular, something similar has been observed'^"^ in the 
hole-doped cuprates. However, the MF theory is problematic, in that the 
Mott gap is associated with long-range Neel order, and the MF model 
predicts anomalously high values for Tjv- 



''Actually, in Ref. a third neighbor hopping t" was included to optimize the fit to the 
experimental Fermi surface curvature. The changes induced by this parameter are small, 
and it will be ignored in the present calculations. 
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2. Mode Coupling Theory 

This anomalous behavior can be cured by incorporating the role of fluc- 
tuations. Indeed, it is known that in a two-dimensional system, the Neel 
transition can only occur at T = - the Mermin- Wagner (MW) theorem^^. 
By treating fluctuations within a mode coupling analysis, the MW theo- 
rem is satisfied^^, and the Mott gap is completely decoupled from long- 
range spin-density wave (SDW) order. Even though = a large Mott 
(pseudo)gap is present even well above room temperature near half filling 

- due to short-range AFM order. The MF gap and transition temperature 
are found to be approximately the pseudogap and T* the onset temperature 
for the pseudogap opening. 

The calculation can be summarized as follows. In a path integral for- 
mulation of the Hubbard model^®, the quartic term is decoupled via a 
Hubbard-Stratonovich transformation and the fermoinic degrees integrated 
out. The resulting action is then expanded to quartic; order in the Hubbard- 
Stratonovich fields 4>. The quadratic interaction reproduces the RPA theory 
of the Hubbard model - the quadratic coefficient is just U5oq, where 6oq is 
the (inverse) Stoner factor Sq = 1 — Uxo{Q + Q,i^) (here Q = (tt, tt) is the 
wave vector associated with the commensurate SDW). The quartic inter- 
action, parametrized by the coefHcient u evaluated aX Q, oj = 0, represents 
coupling between different magnetic modes. The effects of this term can- 
not be treated in perturbation theory, and a self-consistent renormalization 
(SCR) scheme^ is introduced to calculate the renormalized Stoner factor 
Sq = Soq + X. A self-consistent equation for 6 is found, which can be solved 
if the band parameters and the interactions U and u are known. 

To simplify the calculation, A is assumed to be independent of q and oj, 
and the Stoner factor is expanded near Q as 

5q{aj)=6 + Aq'^ -Buj^ -iCuj, (1) 

similar to the form assumed in nearly antiferromagnetic Fermi liquid 
(NAFL)^^ theory. The imaginary term linear in frequency is due to the 
presence of low energy magnon excitations in the vicinity of the 'hot spots' 

- the points where the Fermi surface intersects the Brillouin zone diagonal. 
Here fluctuations toward long-range Neel order lead to strong, Bragg-like 
scattering which ultimately leads to the magnetic Brillouin zone boundary 
at Tjv. Independent of the parameter values, it is found that 5 > for 
T > the MW theorem is satisfied, while for electron doping (5 ^ as 
T ^ up to a critical doping - there is a QCP associated with T = SDW 
order. For hole doping the calculation breaks down - the parameter A is 
negative in a significant doping regime, as will be discussed further below. 
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The key insights of mode couphng theory are: 

• The Mott transition is dominated by hot spot physics, which creates 
zone-edge magnons. The condensation of these magnons creates a new zone 
boundary, and opens up a gap (the Mott gap) in the electronic spectrum. 
In two dimensions (2D), there can be no Bose condensation at finite tem- 
peratures, but the pileup of lowest energy magnons as T decreases leads to 
the appcarence of a Mott pseudogap and a T=0 transition to long-range 
SDW order. 

• Evidence for the existence of local magnons comes from well-defined 

plateaus in the spin susceptibility. Plateaus arc scon in (a) the doping 
dependence of the susceptibility at Q, (b) the q dependence of the suscep- 
tibility near Q, and (c) the uj dependence of the susceptibility (both real 
and imaginary parts) at Q. These plateaus introduce cutoffs in the q and uj 
dependence of the Stoncr factor, Eq. 1, which in general cannot be sent to 
QO, in contrast to the NAFL model. Also, the fiatness of the plateau tops 
makes it difficult to estimate the model parameters from first principles. 
In particular A is strongly temperature dependent, while u is anomalously 
small (this problem had been noted previously^*). 

• A finite Neel temperature can be generated by interlayer coupling. 
In the cuprates, such coupling is typically frustrated, and Neel order more 
likely arises from a Kosterlitz-Thouless transition, after the spin dimension- 
ality is reduced by, e.g., spin-orbit coupling^^. 

3. Results 

3.1. Susceptibility 

In analyzing the ARPES data, a tight-binding band is assumed. 



with a = cos k,a, t = 0.326ey, and t'/t = -0.276. The Hubbard U is 
doping dependent, U/t = 6, 5, 3.5, and 2.9 at a; = 0, -0.04, -0.10, and -0.15, 
respectively, and the mode coupling constant is adjusted to reproduce the 
low-temperature spin stiffness at half filling^", = 0.256ey. The bare 
susceptibility 



efc = -2t{cx + Cy) - At'cxC, 



(2) 





is evaluated near Q to determine the parameters of Eq. 1. 

The susceptibility Xo(Q) 0) has approximately the shape of a plateau as 
a function of doping, Figure la, bounded by the critical points xr and xc 
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Figure 1. (a) Susceptibility xo at Q as a function of doping for several temperatures. 
From highest to lowest curves near x = 0.1, the temperatures are T = 1, 100, 300, 600, 
1000, 2000, and 4000 K. Dashed line = l/U^ff- (b) Circles = pseudo-VHS (peak of xo) 
as a function of temperature Ty; triangles = Tincomm- 




Figure 2. (a) Mean-field transition temperature T* (long dashed line) compared with 
Neel temperature of NCCO and LSCO (solid line) plus magnetic transition temperature 
of stripe ordered phase in Nd-substituted LSCO^^, Circles = model calculation of Tjy 
assuming unfrustrated interlayer hopping, (b) Comparison of mean-field transition (long 
dashed line) to various estimates of pseudogap temperature: from photoemission'^^ (solid 
line), heat capacity^^ (squares) and tunneling^* (circles = A/3, with A the tunneling 
gap). Short dashed line = superconducting Tc- 

where the Fermi surface ceases to have hot spots. These special points act 
as natural phase boundaries for antiferromagnetism, due to the sharp falloff 
in X off of the plateau. 

The special point xh coincides at T = with the Van Hove singularity 
(VHS) of the band. Remarkably, the susceptibility peak has a strong tem- 
perature dependence^^, defining a pseudo-VHS; Fig. lb shows the tempera- 
ture Tv{x), at which the susceptibility peaks at x. This can be understood 
by noting that the energy denominator of x, Eq. 3, is independent oft' , and 
thus would lead to a large peak at half filling, x=0 (associated with states 
along the zone diagonal). At low temperatures, the difference in Fermi 
functions in the numerator cuts this off, and forces the peak to coincide 



6 



1.230 




(1,0.8) (1,1) (0.9,0.9) -0-1 00 xO.1 0.2 



Figure 3. (a) Expanded view of susceptibility xo on the plateaus near Q for a variety 
of dopings at T = IK. From highest to lowest curves near Q, the chemical potentials 
are = -0.25, -0.22, -0.21, -0.20 (x=0), -0.15, -0.10, -0.05, and -0.02eV. All curves 
except /I = — 0.20ey have been shifted vertically to fit within the expanded frame, (b) 
Calculated A{x). 

with the VHS. As T increases, more states near the zone diagonals become 
available, causing the peak susceptibility to shift towards half filling. For a 
temperature Tincomm slightly above Ty, the curvature A becomes negative, 
signalling the instability of the commensurate AFM state. 

The dashed line in Fig. la represents l/Uefj, where C/e// is a doping 
dependent Hubbard U, estimated from a screening calculation^^. The inter- 
section of the dashed line with one of the solid lines defines the mean-field 
Neel transition, XoUeff — 1, Fig. 2a (long-dashed line). As will be shown 
below, once fiuctuations are included, the mean field transition turns into a 
pseudogap onset T* , while the actual onset of long range magnetic order is 
suppressed to much lower temperatures. From Fig. 2b, it can be seen that 
the mean-field T* is consistent with a number of estimates^^'^'^'^"' of the 
experimental pseudogap for hole doping, while a simple calculation of the 
three-dimensional Neel transition associated with interlayer coupling (cir- 
cles in Fig. 2a - see the Appendix for details) can approximately reproduce 
the experimental observations (solid lines) - if the transitions associated 
with magnetic order on quasistatic stripes^^ are included (squares). 

In addition to plateaus in doping, the hot spots lead to plateaus in the 
frequency and wave number dependence of xoiQT^)- For instance. Fig. 3a 
shows plateaus in xoiQ + 0) at a series of dopings at T IK. Once xo is 
known the parameters A and C of Eq. 1 can be calculated; B is quite small, 
and can in general be neglected. The plateau width Qc —> as x ^ xc, 
leading to a strong T-dependence of A, Fig. 3b. 

+ 
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Given the model parameters, the self-consistent equation for S 
becomes^^ 



^2CT^(5 + Ag'2)2+a;2- 
which has the approximate solution 



(4) 



A 3^a2T , 2Cr 



where Sq = 6o + r] — 1, 

3ug2a2,l _„ fan~^(ao)s 

= Aq^/ai^, and (/c and arc wavcnumbcr and (normalized) frequency 
cutoffs, respectively, and the exact form of Z is not required. From the 
logarithm in Eq. 5, 6 must be greater than zero for all T > 0, so there is 
no finite temperature phase transition. At low temperatures, the S on the 
left hand side of Eq. 5 can be neglected, leading to a correlation length 

e = j=eoe''''''/^, (7) 

with S,o = eA/2TC and 

nA\So\ 



Here, 5q is the quantum corrected Stoner factor, which controls the T = 
QCP: there is long-range Neel order at T = whenever (5o < 0, or J/xo > '7- 

3.2. ARPES Data 

From the susceptibility, the contribution to the electronic self energy due 
to one magnon scattering can be calculated. The imaginary part of the 
susceptibility can be written 

where the coupling is approximately g'^xa — 3t//2. Since e is peaked near 
zero when q c:^ Q, ImT, is approximately a (5-function at Approxi- 
mating ImT, = —'7rA'^S{ui — C^+q); then 

A' = ^{5 -do). (10) 
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Figure 4. (a)Spectral function for x = 0.0, k = (tt/2, 7r/2), at T = 100, 500, 1000, 2000, 
3000, 4000, and 5000K (larger splittings corresponding to lower T's). (b) A(T); solid 
line: x = 0.0, dotted lines: x = —0.04, short-dashed lines: x = —0.10, long-dashed lines: 
X = -0.15. (c)/m(E) for x = 0, k = (7r,0) at T = 100 (solid line), 500 (long dashed 
line), 1000 (short dashed line), and 2000K (dotted line). 

The importance of this result can be seen by noting that, by Kramers- 
Kronig, 

i?eE(fc,cj) = , (11) 

- ^Uq 

so that 

G-\k,u;)=oj-^j:-Re^^{k,u;) = - "'^ /"•+^^ . (12) 

The zeroes of G^^ arc identical to the mean field results for long-range 
AFM order^^'^^, except that the long-range gap A — U < niz > is replaced 
by the short-range gap A ~ U^y < >. 

Figure 4a shows the spectral function A{k, uj) — Im{G{k, <.^))/vr for x = 
at fc = {tt/2, 7r/2) at a series of temperatures. The spectrum is split into 
upper and lower Hubbard bands, with a gap approximately 2A. The short- 
range order gap A is plotted in Fig. 4b; it is finite for all temperatures, but 
increases significantly close to the mean-field Neel temperature (arrows). 
The net dispersion for two dopings, x = 0, -0.15, is shown in Fig. 5; it is 
in good agreement with the experimental'' and mean-field^^ results. The 
build up of hot spot magnons is reflected in the growth of Im{Y,{k, u))) near 
Lo = £,g^Q, Fig. 4c (note the logarithmic scale). 

4. Implications for Hole Doping 

One expects, and observes, a certain degree of symmetry between electron 
and hole doping: there is a susceptibility plateau associated with hot spots, 
Fig. 1, which terminates near optimal doping; the termination of magnetic 
order in electron-doped cuprates at a QCP near optimal doping is matched 
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Figure 5. Dispersion relations for (a) x=0 and (b) x=-0.15, at T=100K. Brillouin zone 
directions are F = (0, 0), X = (7r,0), S = {tt,tt). 



in hole-doped cuprates by the observation of a QCP near optimal doping 
which appears to be associated with loss of magnetic correlations^'^; and 
in both cases an optimal, probably d-wave superconductivity is found near 
the QCP. Also, it is expected that U would have a similar decrease with 
doping for either electron or hole doping (see the model calculation, dashed 
line in Fig. 1). 

On the other hand, there are also significant differences, not the least of 
which is the magnitude of the superconducting Tc- Most significantly, there 
is considerable evidence for the appearence of nanoscale phase separation 
- either in the form of stripes^^ or blobs^* - for hole doped cuprates, while 
the evidence is weaker or absent for electron doping^^''^". As noted above, 
this difference arises naturally within the present calculations, which find 
instability of the commensurate magnetic order for hole doping only (e.g.. 
Fig. 3). 

Given these similarities and differences, the present calculations can 
shed some light on aspects of hole-doped physics: 

(1) A large pseudogap is present at half filling, associated with short- 
range magnetic order. Just as for electron doping, it should persist with hole 
doping until short-range magnetic order is lost. The observed^'^ connection 
of the loss of magnetic fiuctuations with the collapse of the pseudogap near 
X = 0.19 strongly suggests an identification of the observed pseudogap at 
T* with the Mott pseudogap, as found for electron doping. While a number 
of theories have proposed that the pseudogap is associated with supercon- 
ducting fluctuations, these seem to turn on at a temperature lower than 
T*.^^. Remarkably, despite the complications associated with stripes, the 
mean-field transition temperature is within a factor of two of the observed 
pseudogap temperatures. Fig. 2b. 

(2) The present calculations point to a close connection between the 
Van Hove singularity (VHS) and the instability of commensurate magnetic 
order. Fig. lb. This strongly suggests that the VHS is responsible for the 
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asymmetry between electron and hole doping, and in particular for any 
frustrated phase separation. 

(3) If the stripes are associated with frustrated phase separation, it is 
important to identify the second (metallic) phase and understand what in- 
teraction stabilizes it (particularly since this is likely to be the phase in 
which high-Tc superconductivity arises). In a purely Hubbard model, this 
would be a ferromagnetic phase, hence probably incompatible with super- 
conductivity. However, the reduction of U with doping found here strength- 
ens the case for a nonmagnetic charge stripe associated with interactions 
beyond the Hubbard model^^. 

5. Discussion: Polaron Limit 

In the very low doping limit, isolated charge carriers should form (spin 
or charge) polaronic states for either sign of doping. The asymmetry be- 
tween hole and electron doping would then be reflected in interpolaronic 
interactions being strongly attractive for hole doping. For electron dop- 
ing, the isolated polarons could be much more easily pinned in the AFM 
background, leading to the much stronger localization found in NCCO^^. 
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Appendix: c-Axis Coupling 

A toy model is introduced to study the effect of interlayer coupling on gen- 
erating a finite Neel transition temperature T^. The interlayer hopping is 
assiimcd to be a constant tz independent of in-plane momentum. While 
a term of the form tz{cx — Cy)^ would not greatly change the results, in 
the physical cuprates alternate Cu02 planes tend to be staggered, which 
should lead to frustration tz{Q) = 0, and greatly reduced interlayer cou- 
pling. Indeed, in the cuprates it is entirely possible that interlayer coupling 
is negligible, and that the Neel transition is actually of Kosterlitz-Thouless 
form, due to reduced spin dimensionality caused by spin-orbit coupling^®. 
Nevertheless, it is instructive to see how constant-tx interlayer coupling 
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can generate a finite T/v. The revised Eq. 5 can be written in the symbohc 
form 

^^ = ^"° + |'"^^)' (^-^^ 

where Tq = ii^ A/Qua? and Dq = 4:CT/e arc (doping-dcpcndcnt) constants 
(Eq. 5 - the extra 7r/2 in Tq coming from the g^-integral) and"^^ D oc 
t^. Thus for finite t^, there is a non-zero Tn given by the solution of 
So + T/T^ = 0, with To* = ro/ln(i:>o/£'). For the calculation in Fig. 2a, a 
constant Tq* = 1200if was assumed, but it is interesting to note that when 
the correct doping dependence of the parameters is included, Tat ^ as 
A ^ Q, suggesting that the much steeper falloff of Tjv with hole doping is 
related to phase separation. 

Equation A.l can be rewritten using Eq. 8. The Neel transition occurs 
when 

JMTN)/UTN)f ^TTn. (A.2) 

with JJJ = {tjtf, r = IQC/edU, d = D/tl and J = At^/U, suggestive 
of a form of inter layer coupling proposed earlier^^. 
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